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Abstract 


The cross-sectional profile of the equatorial electrojet is deter- 
mined for various longitudes along the dip equator, using a set of 4® 
Gauss coefficients for the earth’s magnetic field. The electrojet pro- 
file is found to vary appreciably with longitude due to the asymmetry 
of the magnetic field. To demonstrate the longitudinal variation of the 
electrojet the profiles for longitudes 80° E (for India) and 280° E (for 
Peru) are shown. The maximum current density in the electrojet for the 
latter longitude is found to be much greater than that for the former 
longitude. Also the width of the electro jet is greater at 280° E than 
at 80° E. The magnetic field produced by the electrojet is computed for 
various longitudes along the dip equator and also for the locations of 
six magnetic observatories near the dip equator. The theoretically ex- 
pected field changes are qualitatively in agreement with the observed 
values for these observatories. The average electron drift velocity in 
the electro jet is estimated; over Peru this velocity reaches a few km/sec 
during a sunspot maximum period at the height of the peak of the electro- 
jet. This value is considerably higher than earlier estimates and is 
well above the critical electron velocity, i.e., the ion thermal velocity, 
above which the two-stream instability can occur. 



Introduction 


The recent rocket explorations of the equatorial electrojet in 
India (Maynard and Cahill, 1965a) and off the coast of Peru (Davis et 
al., 1965} Maynard and Cahill, 1965b) have stimulated general interest 
in this subject, and it seems timely to present a model electrojet with 
which the new measurements may be compared. The magnetic observations 
made in the equatorial regions at ground level have been discussed 
extensively (Onwumechilli, 1959a, 1959b; Osborne, 1962a, 1962b, 1963} 
Forbush and Casaverde, 1961). The model presented here is based on a 
model atmosphere described below and on the magnetic field configuration 
determined with the 48 Gauss coefficients given by Jensen and Cain (1962). 
These coefficients were chosen as being the most accurate available at 
the time this work was begun. A recalculation made with a more recent 
model (Hendricks and Cain, 1966) has shown that changes in the final 
results are sufficiently small that the conclusions drawn in this paper 
are not altered. 

Assuming that the physical quantities of the ionosphere are the 
same everywhere along the dip equator except for the differences arising 
from the asymmetry of the earth 1 s permanent magnetic field, the cross- 
sectional profile of the electrojet is determined for different longi- 
tudes. In particular, two such profiles, one for longitude 80°E for 
India and another for longitude 280°E for Peru, are presented. Using 
such profiles the magnitude of the magnetic field due to the electrojet 
as observed at ground level is computed as a function of longitude. The 
longitude variation so obtained is shown to be in fair agreement with 
the observations at several equatorial observatories. 
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The average electron drift velocity in the electro jet is estimated, 
and is shown to be well above the threshold velocity for the two-stream 
instability discussed by Farley (1963a, 1963b) and Buneman (1963). 


Conductivity profile 

We take a rectangular coordinate system (x, y, z) such that: the 


x-y plane coincides with the ground with the positive z axis pointing 


upward} the y axis is along the dip equator its positive direction being 


eastward; the x axis is directed south. The primary electric field is 
assumed to be static and parallel to the y axis; we denote this field 
by Ey. If we further assume that the vertical Hall current is entirely 
inhibited by the (vertical) polarization of the medium, this polarization 
field E z is given by 

£ z = [°2/(°o tan 2 I + a 1 )] Ey 

where a 0 , C]_, and 02 are the direct, the Pedersen, and the Hall con- 
ductivity, respectively, and where I is the magnetic dip (Hirono, 1952; 
Baker and Martyn, 1953). Thus at the dip equator, where 1=0, the 
vertical polarization field is o 2 /cr-]_ times the original horizontal field. 
In Figure 1 the ratio a 2 /oi is plotted against height. (The ionospheric 
parameters used for the calculation of the conductivity elements are given 
in the Appendix.) In the region near 100 km level CT 2 / a l ex tremely 
large, being 20 to 30. It is essentially this large vertical polarization 
field that dirves the electrojet. The sum of the Pedersen current Ey 

and the Hall current o^E z is the total electrojet current. At the dip 

2 

equator this is o^Ey, where 4- c 2 /o^; is called the Cowling 

conductivity. 
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We denote the effective conductivity in the y direction by a , 

where 

2 2 

a = °1 + °2 / (c o tan 1 + a i) 

Since the electric field E is assumed to be uniform, the cross-sectional 
profile of Gyy in the x-z plane gives the corresponding profile of the 
relative current density. The ionospheric parameters are taken to be 
functions of z alone; while the magnetic field as determined from the 
48 Gauss coefficients varies with longitude along the dip equator. It 
is only the latter that generates the longitudinal variation of the 
electrojet profile. The electron density distribution at noon is not 
likely to be the same everywhere along the dip equator, but what we in- 
tend to show here is the effect of the nonuniformity of the magnetic 
field along the dip equator, separating this from other possible effects. 
We make an idealization that the primary electrostatic field is assumed 
to be uniform and parallel to the dip equator. The exact electric field 
distribution can be determined only if the wind pattern and the distri- 
butions of the electron density and other ionospheric quantities are 
completely known all over the world. Without a complete knowledge of 
these factors we have to rely on known properties of the magnetic vari- 
ation Sq and of the associated electric field theoretically deduced. 

The general direction of the electric field in midday in the equatorial 
region as inferred by the dynamo theory is approximately eastward (Maeda, 
1955) . The characteristics of Sq are known to be similar along the lines 
of equal dip, and near noon the Sq current in low latitudes appears to 
flow nearly parallel to the dip equator (Price and Wilkins, 1963). Thus 
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for an idealized model the assumption that the primary electric field 
is parallel to the dip equator seems to be appropriate. 

Figures 2 and 3 show contours of equal values of the effective 
conductivity Oyy in the vertical plane normal to the dip equator for 
longitude 80°E (for India) and for longitude 280°E (for Peru). The 
numbers indicated along the contours are in units of 10”^-5 e.m.u., or 
10”4 mho/meter in m.k.s. units. With the assumption that Ey is uniform 
the electric current density is proportional to <Jyy, hence the profiles 
can be thought of as the current density profiles. It should be noted 
that the horizontal scale is 20-fold contracted) therefore, the actual 
profiles are much more oblong horizontally than is seen in these figures. 
Electro.jet over Peru (280°E) 

The peak of Cyy is at 103.5 km altitude directly above the dip 
equator, and its value is 118 in the units used in Figure 3> i.e., 

1.18 x 10 -W e.m.u. (1.18 x 10~^ mho/meter). The altitudes at which Oyy 
is reduced to one-half the peak are 97 km and 110.5 km directly above 
the dip equator. The contour for this half -peak value extends horizontally 
from 50 km north to 80 km south of the dip equator. It is noted that 
the electrojet profile is not symmetric with respect to the dip equator) 
this feature is solely due to the asymmetry of the magnetic field with 
respect to the dip equator. 

Electro.jet over India (80°E) 

The Oyy peak is at 102.0 km altitude above the dip equator, and its 
value is 71 in the same units as in Figure 2, i.e., 0.71 x 10 -1 3 e .m.u. 
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(0.71 x 10~ 2 mho/meter). Thus over India the peak altitude is about 

1.5 km lower than over Peru and the peak value over India is only 60 % 
of that over Peru. The altitudes at which Oy y is one-half the peak are 

95.5 km and 110 km. Thus the thickness of the electro jet measured by 
these altitudes is nearly the same in India as in Peru; the thickness so 
defined is about 14 km. However, the steepness of the conductivity 
variation with height is much more pronounced over Peru than over India. 
Horizontally the contour for the half-peak value extends from 55 km 
north to 40 km south of the dip equator over India. There is a tendency 
that the electrojet stretches more to the south than to the north in 
Peru, while in India the tendency is opposite. It is noted that not 
only is the peak current density smaller over India than over Peru the 
width of the electrojet is much narrower over India. To illustrate the 
large differences between the profiles the contours for these longitudes 
are superimposed upon each other in Figure 4. 

It is interesting to note that the altitude of the major part of 
the electrojet is a little lower than the altitude at which the normal 
Sq current flows at latitudes slightly off the dip equator as if the 
electro jet hangs below the Sq current sheet. 

At other longitudes 

As has been shown above for two representative longitudes the peak 
value of CTyy over the dip equator varies considerably with longitude. 

The ma ximum of this Oyy peak occurs near longitude 310°E where its value 
is 1.33 x 1CT 13 e.m.u. and the minimum near 90°E where the peak value of 
is 0.68 x lO - -'- 3 e.m.u. Thus with regard to the Oyy. peak India is 
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close to the minimum and Peru to the maximum, making this pair well 
suited for further experimental studies of longitudinal inequality of 
the electrojet. 

Having constructed conductivity profiles for different longitudes 
it is possible to calculate the longitudinal variation of the magnetic 
field produced by the model electrojet. The horizontal magnetic field 
AH 0 at the dip equator due to the electrojet is shown in Figure 5 as a 
function of longitude} values are normalized in the manner explained 
below to make them comparable with the observational data. The scale 
can be adjusted, however, by multiplying a constant. It is noted that 
the value of AH Q at 280°E is almost twice that at 80°E. 

To compare the computed field with the observations at the six 
magnetic observatories listed in Table 1 the horizontal component of 
the electrojet field at the actual locations of these observatories was 
computed. There will be an induction effect from the conducting earth, 
and this effect may be different at different observatories. However, 
for the sake of simplicity we assume that the horizontal component of 
the magnetic field due to the electro jet is increased by the induction 
at the same rate everywhere. Then this factor can be absorbed in the 
normalization constant in the computed field. For the observational 
data we take the peak amplitude in H measured from the midnight value 
for the annual average H variation for the five international quiet 
days for 1958. Denoting these values by the observed AH, the normalization 
constant for the calculated AH is determined so that the average of the 
latter over the six observatories becomes equal to the corresponding 
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average of the observed AH. The observed and the calculated Ah so 
defined are plotted in Figure 6. 

The agreement between the computed and the observed values of AH 
is only qualitative. The computed Ah for Huancayo is too large and 
those of Trivandrum and Ibadan too small compared with the corresponding 
observed values. The longitudinal variation of AH calculated here seems 
to be much larger in amplitude than is observed. Without quantitative 
arguments we point out that the change in the conductivity with long- 
itude will reduce the longitudinal variation of the electrojet intensity 
because of the electrostatic field set up by the electrojet itself; in 
particular, at longitudes near the conductivity maximum at 310°E long- 
itude the electrojet may be appreciably suppressed by this effect. It 
is also pointed out that the observed values used in this paper are the 
annual averages of the amplitude in H and are only crude measures. The 
primary purpose of the present paper is to present a model profile of 
the electrojet, and not to discuss the longitudinal variation in the 
intensity of the electrojet in precise terms. Also the observed average 
amplitude Ah (measured from the midnight level) has large seasonal vari- 
ations. The amplitude is largest in equinox for all the stations except 
for Jarvis Island where it is a maximum in December solstice; while the 
amplitude is smallest in December solstice for Ibadan, Addis Ababa, and 
Trivandrum, but for Jarvis Island, Koror, and Huancayo the minimum is 
in June solstice. Thus it appears that factors other than the longi- 
tudinal inequality in the permanent magnetic field play an important 
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part in producing the longitudinal variation in the electrojet. 

As to the width of the electrojet there is little observational 
data to compare with the present result. However, Onwume chilli (1959b) 
estimated the width of the electrojet in Nigeria to be about 440 kmj 
while Forbush and Casaverde (1961) took the electrojet width to be about 
660 km in Peru. Since these figures are based on different models, they 
cannot be compared directly. Nevertheless the tendency that the elec- 
trojet is of greater width in South America than in Africa is in complete 
agreement with our present results. 

Magnetic field distribution through the electro.let 

To provide theoretical values to be compared with the recent rocket 
measurements the magnetic field was computed as a function of height 
directly above the dip equator both for longitudes 80°E and 280°Ej the 
results are shown in Figure 7. In the calculation the value of Ey was 
adjusted so as to give values of AH at ground level comparable with the 
typical observed values in India and in Peru during the periods the 
rocket firings were made. For 80°E, i.e., for India, the electrojet 
field at ground level on the dip equator is 24VJ if "the induction effect 
is taken to be one-half of this value, the total field will be 36v. For 
280°E, i.e., for Peru, the electrojet field at the ground is taken to be 
53 y s° that with the same rate for the induction effect as above the 
total field would be about 80y. These values are very much lower than 
those observed in 1958 because the solar activity is near a minimum 
during the period we are concerned here. In Figure 7 the magnetic 
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field directly due to the electrojet alone is shown without the induced 
field. When the induction effect is included the magnetic field gra- 
dient below 90 km would tend to become smaller than is shown in these 
figures. 

With the present model there is a sudden reversal in the vertical 
gradient of the magnetic field at about 90 km altitude, above which the 
field change is very nearly linear to about 115 km. There is a second 
turning point at about 125 km to 130 km altitude, and the change becomes 
nearly linear again above this level. 

Discussions 

1. Rastogi (1962) compared the observed electrojet effects in the 
American, African, and Indian zones and showed that the effect is most 
pronounced in the American zone, less so in Africa and least in India. 
This is in agreement with the present result shown in Figure 5. By an 
approximate calculation Rastogi argued that at 100 km level varies 
very little with longitude, and suggested that the observed large long- 
itudinal variation is probably due to a longitudinal variation in the 
electric field. However, the present calculation shows that varies 
appreciably with longitude, the ratio of at 80°E to that at 280°E 
being 0.61 at 102.5 km altitude. The errors introduced by the approxi- 
mations in Rastogi' s argument contribute enough to mask the dependence 

2 , 

of 0£ on the magnetic field. The wind pattern and the electron 
density distribution are not likely to be exactly the same along the 
equator, but we point out that the asymmetry of the magnetic field alone 
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could give rise to a considerable longitudinal variation in the elec- 
trojet intensity. 

2. In the dynamo theory the conductivity is usually integrated 
over the thickness of the ionosphere. The integrated value of our a ^ 
from 90 km to 160 km altitude is compared in Table 2 with the height- 
integrated a 3 given by other workers. Although these values are in 
agreement in magnitude, it seems that for detailed studies of the 
global Sq current system the variation in the conductivity with long- 
itude (as well as with latitude), in addition to the variations with 
season and solar cycle, should be taken into account. A new calculation 
is under way to construct an electro jet model based on the actual rocket 
measurement of electron density recently made over the dip equator off 
the coast of Peru. 

3. Since the Hall current due to the vertical polarization field 
is the main constituent of the electrojet, the average electron velocity, 
v g , can be approximately represented by 

v e = (E z /B)/(1 + v//^) 

where v e and u) e are the electron collision and cyclotron frequencies 
respectively. The direction of the electron stream is westward. At the 
electro jet height v e « u) e , and hence the average electron velocity is 
very nearly E z /B. Therefore, if Ey is uniform, v Q varies with height 
proportionately with the ratio a 2 /<?]_• Thus when Ey/B is multiplied to 
the abscissa in Figure 1, it approximately gives the average drift 
velocity of the electrons in the electrojet els a function of height. 
Throughout the height range 90 km to 130 km over Peru v e is given, to a 
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good approximation, by 3.8 x 10^- Ey ®/sec when E^. is expressed 

in volts/m. When E = 2.4 x 10“3 volts/m, AH at Huancayo due directly 

J 

to the electrojet is lOOyj if the induction effect is taken to be one- 
half of this external field, the total variation observed will be about 
150\. With this value of Ey, the average electron drift velocity at 
the height of the peak current is 2.7 km/sec. The average ion velocity 
is little affected by the presence of the large polarization field in 
the electrojet, and is only of the order of a few meters per second at 
the most in the 100 km to 110 km region. During sunspot minimum years 
the a-verage electron streaming velocity at the electrojet peak would be 
about one-half of the above estimate. 

Farley (1963a, 1963b) and Buneman (1963) have shown that if the 
streaming electron velocity relative to the corresponding ion velocity 
exceeds the ion thermal velocity a two-stream plasma instability can 
occur and they attributed a certain type of E region irregularities 
observed at the dip equator (and in the auroral zone) to longitudinal 
ion waves generated by this instability. The observational aspect of 
this interesting phenomenon has been discussed in detail by Bowles, 
Balsley, and Cohen (1963) and Cohen and Bowles (1963) . The ion thermal 
velocity in the electrojet region is about 400 m/sec. If we take this 
value as the minimum critical electron drift velocity for the instability, 
the present model over Peru shows that if Ey is taken to be 2.4 x 10 J 
volts/m, the average electron drift velocity in the electrojet region is 
well above this critical value. 
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Zmuda (i960) estimated the electron drift velocity at 310 m/sec. 

This value is based on a current density of 1 x 10“5 amp/m 2 with an 
electron density of 2 x 10^ electrons/cm^. His estimate of the average 
electron velocity is considerably lower than our estimate. At the 
height of the peak current density in the present model the electron 
density is 8.0 x lcA electrons/cm^, and the current density is 
2.8 x 10“5 amp/m 2 . As was pointed out by himself the high electron 
density Zmuda used corresponds to the density in a sporadic E layer; 
while the present electron density may be thought to be a value under 
a normal condition. In Zmuda' s model the electric current density is 
averaged over the thickness of the electrojet; therefore, the current 
density quoted above must be considerably lower than that at the peak 
of the electrojet. As to the estimate of the electric field, Zmuda' s 
estimate, 4 x 10“4- volts/m, is lower than the value adopted here, namely, 
2.4 x 10"3 volts/m. This is because Zmuda took a ^ to be constant through 
the electroj et legion ; moreover his value of 2.5 x 10 mho/m, 

is higher than the present estimate at the electrojet peak, i.e., 

1.18 x 10“ 2 mho/mj his is high because of the high electron density 
in his model. Matsushita (1962) estimated the horizontal electric field 
to be 1.6 x 1CT 3 volts/m, and obtained an estimate of the average electron 
drift velocity of 500 m/sec by averaging the current over the thickness 
of the electrojet. Since these authors assumed the current to be uni- 
form over a thickness of a few tens km, the current density is greatly 
reduced compared with the present current density estimate for the peak 
electroj et altitude. 
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APPENDIX 


The conductivity elements o Q , oj, and 02 (in e.m.u.) are calculated 
by the following formulas: 

c 0 = N e e 2 [(mgVg) -1 + (m^)- 1 ] 

a l = N e e2 Um e v e ) _1 v e 2 (v e 2 + u>e 2 ) -1 + (m-jVi) -1 + uui 2 )" 1 ] 

a 2 = N e e 2 [(m^)- 1 v e co e (v e 2 + u> e 2 ) -1 - (v ± 2 + (Dj 2 ) -1 ] 

where 

N e = electron density, electrons/cm^ 
e = magnitude of the electronic charge 
m 6) i = electron, or ion, mass 
v e ^ = electron, or ion, collision frequency 

= electron, or ion, cyclotron frequency, taken to be positive 
regardless of the sign of charge 

It is assumed that there is only one type of singly ionized ions. 

Because of the polarization of the medium the electron density is 

N 

not exactly the same as the ion density. However, the deviation (AN e ) 
from the neutrality is extremely small. Assuming that all the quantities 
concerned vary with z alone, we can estimate the charge density from 
Poisson's equation: dE z /dz = ^hc^eANg. 

Thus 

AN e = Uttc^)"" 1 Ey d{c 2 /o 1 )/ dZ 

Taking Ey = 2.4 x 103 e.m.u., AN e is of the order of 10 “4 cm -3 in t }, e e i ec _ 
trojet region. Thus the assumption of neutrality in the computation of 
the conductivity elements is justified. 
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Following Chapman (1956) the collision frequencies are computed by 
the following formulas: 

v e = v ei + v en 

where 

v ei = [34 + 4.18 log (T 3 N e _1 )] N e T -3 / 2 

v en = 5.4 x io-io NnT i/2 

and 

v i = v in = 2.6 x 10-9 ( Ne + Nn ) w-!/ 2 
where N n , T, and W are the neutral particle density, the temperature, 
and the mean molecular weight of the ions (which is assumed to be equal 
to that of the neutral particles) . 

The atmospheric model given by Champion and Minzner (1963) was used 
for N n , T, and W. These data together with the computed collision 
frequencies are given in Table 3 in an abbreviated form. The profile of 
the electron density N g used is shown by the heavy curve in Figure 8} 
this curve is a smoothed average of the two dotted curves that represent 
the rocket data obtained at Wallops Island by Bauer and Jackson (1963) 
with the CW propagation technique. The electron density profile adopted 
here is meant only to serve as a model. In a new model to be presented 
in a later publication the electron density determined by the actual 
rocket measurements in the electro jet region and a more recently deter- 
mined spherical harmonic representation of the main field will be used. 
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Table 1. Six equatorial observatories used in this paper. 


Distance from dip 


Observatory 

Longitude 

Latitude 

Dip* 

equator to station** 

Ibadan 

3 . °90E 

7.°43N 

-5.°29 

233 km 

Addis Ababa 

38.77E 

9.03N 

-1.42 

56 

Trivandrum 

76.95E 

8. SON 

-1.94 

83 

Koror 

134. 50E 

7.33N 

-0.22 

21 

Jarvis Island 

160. 03W 

0.38S 

1.77 

-80 

Huancayo 

75.33W 

12.05S 

0.70 

-48 


* Computed values based on the spherical harmonic representation of the 
main field, and not observed values. 


** The distance between the computed dip equator and the station; positive 
when the station is south of the dip equator. 



Table 2. Comparison of estimates 
several workers. 

present paper 80°E 
280°E 

H. Maeda (1953) 

Baker and Martyn (1953) 

Fejer (1953) 

Chapman (1956) hot ionosphere 
cold ionosphere 
Hirono and Kitamura (1956) 


of the height-integrated by 

^ 2i x lO - ^ e.m.u. 

1.97 

0.775 

1.64 

1.80 

2.08 

1.93 

1.52 


Table 3. The atmospheric model used in this paper 


Height 

N n 

W 

T 

v e 

v i 


cm“3 


°K 


sec 

1 

S< 

90 ]m 

6.584 x 10 13 

29.00 

h 

180 

4.77 

x 10 5 

3.18 x 10 4 

100 

1.036x 10!3 

28.88 

210 

7.61 

x 10 4 

5.02 x 103 

110 

2.073 x 10 12 

28.56 

257 

1.63 

x 10 4 

1.01 x 103 

120 

5.235 x 10 11 

28.07 

349 

4.72 

x 103 

2.57 x 10 2 

130 

1.650 x 10 11 

27.70 

515 

1.85 

x 103 

8.15 x 10 1 

140 

7.476 x 10 10 

27.31 

715 

1.03 

x 10 3 

3.72 x 10 1 

150 

4.115 x 10 10 

26.93 

893 

7.22 

x 10 2 

2.06 x 10 1 

160 

2.622 x 10 10 

26.65 

1022 

5.89 

x 10 2 

1.32 x 10 1 



Figures 


Figure 1. 
Figure 2. 


Figure 3. 
Figure 4* 
Figure 5. 
Figure 6. 
Figure 7. 
Figure 8. 


The ratio (f° r 280°E longitude) as a function of height. 

The conductivity (°yy) profile for longitude 80°E for India 
in the vertical plane perpendicular to the dip equator; con- 
tours of equal a are drawn in units of 10 e.m.u., or 

tJ «/ 

10 -/ + mho/meter; north to the right. Note that the horizontal 
scale is 20-fold contracted in comparison with the vertical 
scale. 

The conductivity ( CT yy) profile for longitude 280°E for Peru, 
units are the same as in Figure 2. 

A comparison of the conductivity profiles for longitudes 80°E 
and 280°E; units are the same as in Figures 2 and 3« 

The computed electrojet field at the dip equator as a function 
of longitude. 

Comparison of the observed values of AH with those computed, 
for six equatorial observatories. 

The magnetic field produced by the electrojet as a function 
of height, for longitudes 80°E and 280°E. 

The electron density profile used in this paper. Two curves 
drawn with broken lines represent profiles obtained by 
rockets over Wallops Island (Bauer and Jackson, 1963). 
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